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Abstract 

The quantum theory of hnearized perturbations of the gravitational field of 
a Schwarzschild black hole is presented. The fundamental operators are seen 
to be the perturbed Weyl scalars and ^4 associated with the Newman- 
Penrose description of the classical theory. Formulae are obtained for the 
expectation values of the modulus squared of these operators in the Boul- 
ware, Unruh and Hartle-Hawking quantum states. Differences between the 
renormalized expectation values of both |^'o| and 1^*41 in the three quantum 
states are evaluated numerically. 



I. INTRODUCTION 



In this paper we shall study quantized, linear perturbations of the gravitational field of 
a Schwarzschild black hole. Since this represents a study of the one-loop approximation to a 
theory of quantum gravity, it should provide useful insights into what the full theory should 
look like. In addition, there is a growing body of evidence which suggests that the influence 
of non-conformally invariant quantum fields, foremost of which is linearized gravity, strongly 
dominates that of conformally invariant quantum fields in the neighborhood of a space-time 
singularity. For example it has been shown recently [|l| that the value of the renormalized 
energy density for the graviton field in the vicinity of a conical singularity in an otherwise 
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flat space-time (the idealized cosmic string) is roughly ten times that for the electromagnetic 
field and one hundred times that for the conformally coupled scalar field. A knowledge of 
how gravitons behave near the singularity at the center of a black hole is therefore also likely 
to be crucial to our understanding of quantum gravity. 

The classical theory of linearized perturbations of a Schwarzschild black hole was de- 
velopped by Regge and Wheeler, 0] Zerilli and others. Here we follow the approach of 
Teukolsky, |^ who gave a complete set of solutions to the perturbation equations for a Kerr 
black hole within the Newman-Penrose formalism. In this formalism the system is described 
by two quantities — the perturbed Weyl scalars and if 4 (here and in the following, a dot 
indicates a linearly perturbed quantity). These are the field variables of interest describing 
the semi-classical theory, since 

(a) they carry in their real and imaginary parts information on all the dynamical degrees of 
freedom of the perturbed field; 

(b) they are invariant under both gauge transformations and infinitesimal tetrad transfor- 
mations which leave the metric perturbation intact, as befits a physical quantity; 

(c) they are simply expressed (albeit in a particular gauge) in terms of the metric perturba- 
tion which provides the most direct route to quantization; and 

(d) they directly measure the energy flux of classical perturbations of the black hole across 
the horizon and at infinity. [^[^ 

In this work our main concern will be to extend the work of Candelas et ai, P| by first 
presenting formulae for the expectation values of |\E'o| and |\E'4| with respect to the three 
physically relevant quantum states for the Schwarzschild black hole, namely the Boulware, 
Unruh and Hartle-Hawking states, and then evaluating these expressions numerically. Since 

I ' 1 2 I ' 1 2 

|\E'o| and I'^^l are quadratic in the field variables their expectation values will be infinite; 

I ' 1 2 I ' 1 2 

however differences between the expectation values of either |^o| or |\E'4| with respect to 
the various quantum states mentioned above are finite and these differences are the objects 
which we will calculate. (Indeed, some would argue that such differences are more likely to 
be of physical significance than the results obtained from renormalization procedures.) In 
the process of achieving these main goals we also fill several gaps in the literature concerning 
quantized perturbations of black holes; in particular we present both low-frequency analytic 
formulae and a set of Wronskian relations for the reflection and transmission coefficients 
associated with the theory. 

The format of the paper is as follows: the two sections immediately following this in- 
troduction summarize the essentials of what is already known about the classical and semi- 
classical theories of linear perturbations of black holes. Thus in Sec. ||we present Teukolsky's 
complete set of classical solutions for the perturbed Weyl scalars. Since Schwarzschild space- 
time is static and spherically symmetric, the temporal and angular dependences of each of 
these solutions are given by exp(— ia;t) and spherical harmonic functions respectively; the 
radial component however is a solution of a complex, second order, ordinary differential 
equation for which no closed-form general solution exists. We discuss the behavior of the 
radial component both as one approaches the event horizon and spatial infinity. We end 
Sec. m by writing down formulae for the classical energy flux due to the perturbations across 
the event horizon and at infinity. 

Sec. |ITl| is a review of the semi-classical theory. We present the complete set of solutions 
for the metric perturbations due to Chrzanowski, H and verify that these yield Teukolsky's 
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complete set of perturbed Weyl scalars. The question of orthonormality of the mode set is 
rigorously addressed, the perturbations are quantized and expressions for differences between 

I ' 1 2 I ' 1 2 

the expectation values of \'^q\ and |\E'4| in the Boulware, Unruh and Hartle-Hawking states 
are presented in terms of the radial functions. 

In Sec. |IV| we prepare the expressions obtained for the various expectation values for 
numerical evaluation. The first step is to derive a power series representation for the general 
solution to the radial equation; the particular radial solutions occuring in the expectation 
values may then be substituted for these power series, weighted by appropriate constants 
called reflection and transmission amplitudes by virtue of the analogy with classical scatter- 
ing. 

Before the differences between the expectation values can be evaluated numerically, the 
reflection and transmission amplitudes must be determined explicitly. This is achieved in 
Sec. |I^ by comparing the power series representations for the particular radial solutions 
at large values of the radial coordinate with their known asymptotic forms. Graphs are 
presented for these amplitudes against frequency and angular number, and the results are 
checked for both internal and external consistency; the former by means of Wronskian 
relations derived between the coefficients in Appendix and analytic formulae for the 
amplitudes valid at low frequencies derived in Appendix ^; the latter by reproducing Page's 
0] results for the luminosity due to graviton emission from the black hole. 

In Sec. we give an analysis of the asymptotics of the expressions for the differences 
between the expectation values of the perturbed Weyl scalars. 

In Sec. [VI I| we present graphs for these differences, and discuss the main difficulties 
encountered in the numerical computation. 

Throughout we use geometrized Planck units {G = c = h = 1) and follow the sign 
conventions of Misner, Thorne and Wheeler. [110 



II. REVIEW OF THE CLASSICAL THEORY 

For space-times possessing a high degree of symmetry, it is often possible to reduce 
considerably the number and complexity of the equations of general relativity by projecting 
them onto a null, complex tetrad which encodes that symmetry. This is the motivation 
which lies behind the Newman- Penrose (NP) formulation of general relativity. [|lT 



In the case of a Schwarzschild black hole of mass M, whose space-time is described by 
the metric 



g,,^ = diag 



' — , --,r ,r sm 9 

A 



A = r(r - 2M) 



in terms of the usual spherical polar co-ordinates {t, r, 6, 0), an appropriate choice of tetrad 
which mirrors the temporal and spherical symmetry of the space-time is the Kinnersley 
tetrad 



\2Y 



1 / A ^ 



(2.1) 
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e(3,'- = e,r = ^(0,0,l,^). 

The simplification afforded by using this tetrad is exemphfied by the fact that all spin 
coefficients vanish apart from 

1 A M ^ cote , , 

P = — P = n 7 = — ;t — a = p = — 1=- (2.2) 

^ r 2r3 ' 2r2 ^ 2V2r ^ ' 

and all tetrad components of the Weyl tensor are zero apart from 

M 

^2 = ^- (2.3) 

(More than one set of conventions exists for the NP formalism; those employed in this paper 
are spelled out in Appendix |^ to avoid any confusion.) 

Teukolsky took full advantage of the simplifications afforded by the NP formalism 
when investigating perturbations of the gravitational field of a Kerr black hole. He obtained 
linearly perturbed versions of the NP analogues of the basic equations of general relativity 
for the black hole metric and, by working in the Kinnersley tetrad, was able to decouple 
the equations for the ingoing and outgoing radiative parts of the perturbed Weyl tensor, 
\E'o and \E'4. These were found to satisfy the "master perturbation equation," which, in the 
Schwarzschild limit, reads 

aW- + A at - ^8^^ - + 1"'' - 



9^ 19^ 2is cos9 d a in 
cot a— t^tt: 77777 ^ 77" + s cot 



de^ de sin' ^902 sin'^ 



$s = (2.4) 



where s = ±2, $2 = ^0 and $_2 = r^\E'4. It is remarkable that this equation can also be 
shown to describe the behavior of a test scalar field (s = 0, $0 = V^) or electromagnetic field 
(s = ±1, $1 = 007 ^-1 = r'^(p2, where 0o and 02 are the tetrad components F^ye(i)'^e(z)^ and 
Ffiue(4)^e(2)'^ respectively of the Maxwell tensor F^,^) on the Schwarzschild background. For 
this reason we shall keep s arbitrary whenever possible in the ensuing discussion. 
By separating variables in (|2.4| ) the following complete set of solutions is obtained: 



^, = e-'^\RUr),Yr{e,(j)) (2.5) 

where u G [0, 00) and /, m are integers satisfying the inequalities / > |s|, —l<m<l. 
sYi"^{6,(j)) is a spin-weighted spherical harmonic, whose relevant properties are reviewed in 
Appendix B. sRiuj{r) satisfies the ordinary differential equation 



dr \ dr / A 



,RUr) = 0. (2.6) 

The solutions of (|2.6|) cannot be obtained analytically in closed form (we shall solve it 
numerically later). However much can be learned about their properties by exploiting the 
analogy of ( |2.6|) to a classical scattering problem, which becomes apparent when we write 
(PD in the form [| 
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2 + sVluj 



sQiw — 0. 



Here r* is the tortoise co-ordinate 



r + 2Mln 



2M 



- 1 



sQiLuii") = Aar sRiui{r) and sVJa;(?") is the (complex) potential 



5 2iscj(r - 3M) s^M^ r - 2M 
sV^.(r) = + '—^ -— ^ 

As r — > oo the potential ( p.8|) can be approximated by 



/(/ + !) + 



2M 



sViJr] 



2iujs 



so that ( |2.7| ) possesses solutions which behave like 

sQUr) ~ r±^e^^-^ 

which gives 

Similarly, in the limit as r — ^ 2M 



or 



IS 



SO that 



and 



4M 



or 



(2.7) 



(2.8) 



(2.9) 



(2.10) 



f2.111 



(2.12) 



To end this section, we justify the assertion made in the Introduction that the per- 
turbed Weyl scalars \l/o and \I/4 measure the energy flux of classical perturbations of the 
Schwarzschild black hole across the horizon and at infinity. 

Hartle and Hawking ^ have shown that the energy flux transmitted across the horizon 
of a fluctuating classical black hole is given by the formula 



2 rphoT 



dtdfi 



M2 

— \ani2M)l 

TT 



(2.13) 



where CTnir) is the perturbed shear of the null congruence ^e(i)^^^ which generates the 
future horizon, the calculation having been performed in the Hartle- Hawking tetrad which 
is defined as follows in terms of the Kinnersley tetrad ( p.l|) : 
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A 
2^ 



6(1)^ 



2r2 



(e(3)0, = (e(4)-), = e(3)^ 



It can also be shown from the perturbed Newman-Penrose equations that for perturbations 



of frequency u for which 



(Th(2M) 



hm 



'A 



1 + AiMio) 



(2.14) 



(where \l/o is with respect to the Kinnersley tetrad). We thus obtain the formula 



lim 

dtdn 2%M4 (16M2tu2 + 1) r^2M 



A^l^ 



Oa;| 



(2.15) 



At spatial infinity the space-time becomes flat and one can therefore carry over standard 
results from the theory of gravity linearized about Minkowski space, enabling one to derive 
the formulae M 



2 T-iin 



d^E 



1 



dtdn 647rcu2 

^2^out I 



dtdn 47rcj2 



lim 
lim 



Owl 



(2.16a) 
(2.16b) 



for the incoming and outgoing energy flux infinitely far from the black hole. Formulae ( p. 151) , 
( ^.16| ) demonstrate that |^i/o|^ aiid {"^^f directly measure the energy flux of gravitational 
perturbations of a Schwarzschild black hole across the horizon and at infinity. This physical 
interpretation, together with the other features already discussed, makes \E'o and \E'4 the 
natural choice as field variables for the semi-classical theory. 



III. REVIEW OF THE SEMI-CLASSICAL THEORY 

Teukolsky has provided us with a complete set of perturbed Weyl scalars \&o and 
(henceforth collectively labelled \1/a)- One could quantize the theory directly in terms of 
this set 0, however we choose to follow Candelas et al by working in terms of the metric 
perturbation. Thus we require a complete set of solutions for the metric perturbations hfj^„. 
The relationship between and h^i, is obtained by perturbing the defining equations for 
the Weyl scalars. Consider for example 

^0 = -C(i)(3)(i)(3) = -C^„pxe(ife(3)''e(ife(3)^ (3.1) 
which when perturbed linearly becomes 

^0 = -Ct.upxe{i)^e(^3fe^ife^3)^ - C^^px (e(i)''e(3)'^e(i)''e(3)^ 

+e(i)'^e(3)"e(i)^e(3)^ + e(i)^e(3)''e(i)''e(3)^ + e(i)^e(3)"e(i)''e(3)^) . (3.2) 

Since the tetrad spans the tangent space we can write e(^a)^ = A(^a)^''^ ^(^h)^ for some functions 
y4((j)(''); it follows that the last four terms on the right hand side of (|3.2| ) vanish, as the 
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only non-zero, unperturbed Weyl scalar is \E'2 = — C'(i)(3)(4)(2). In the first term one can use 
expressions provided by Barth and Christensen [jl3l for the perturbed Riemann and Ricci 
tensors to express C^^px in terms of h^^i, and (after some tedious algebraic manipulation) 
arrive at the following formula for in terms of h(^a){b)- 

^0 = 1 {{S + 2a)5/i(i)(i) - {2D - 3p){5 + 2a)/i(i)(3) + (D - p){D ~ p)/i(3)(3)} • (3.3a) 
A similar argument for \1'4 yields 

^4 = ^ {{S* + 2a)(5*/i(2)(2) - (2A + 3/i + 47)(<5* + 2a)/i(2)(4) + (A + /i + 27)(A + /i)/i(4)(4)} • 

(3.3b) 



Despite appearances Eqs. ( |3.3| ) can be inverted to yield /i^^, in terms of following a 
procedure due to Wald |]14[ which is facilitated by a choice of gauge reflecting the symmetry 
of the background space-time. It is then straightforward to obtain a complete set of solu- 
tions to the equation of motion for /i^i, from Teukolsky's complete set of Weyl scalars "^a- 
Rather than repeat this lengthy derivation here however, we simply state the result (origi- 
nally obtained by Chrzanowski P] using less direct methods) and verify that it does indeed 
constitute a complete set of solutions to the perturbed Einstein equations, by substituting 
into ( |3.3| ) and reproducing Teukolsky's solution set (pl5|). 

We write Chrzanowski's complete, complex mode set as 

{h^^il, m, u, P- x), h^:{l, m, u, P; x)}^_^^^_^_^ (3.4) 
where A G {in, up} and P = ±1. The explicit form of the modes for which A = in is 

h}^M,m,u;,P;x) = AT- {e^. +2l^"^(e, 0) + PB;, _2l^™(e, 0)} -2RZir)e-'^' (3.5a) 
in the ingoing radiation gauge h^^yCii)' = 0, h^'^ = 0; when A = up we have 

h^lil, m, P; x) = N^^ {T,, ^^Yne, <P) + PT;, +2l^'"(^, <P)] +2PrJ(r)e-'-* (3.5b) 

in the outgoing radiation gauge h^^e(^2)'^ = 0, h^,'^ = 0. Here A^'^ are constants which will 
be fixed by the quantization prescription, and Qp,u, '^nu are the second-order differential 
operators: 

= -e(l)^,e^l)„{S* - 2a) (5* - 4a) - 6(4)^6(4)^ (P* - p)(P + 3p) 
+ ^ (e{i)/.e(4)^ + 6(4)^6(1)^) [D{6* - 4a) + {6* - 4a) (P + 3p)] 
'^tiu = P~'^ {-£(2)^*6(2),. (5 - 2a){6 - 4a) - e(3)^e(3)^(A + 5p - 27)(A + p - 47) 

+ \ (e(2)/.e(3). + 6(3)^,6(2).) [{6 - 4a)(A + /i - 47) + (A + 4/i - 47) (5 - 4a)] | . (3.6) 

The quantities _2P|Ji(r) and +2P/^(^) appearing in ( p.5| ) are those particular solutions 
of the radial equation (|2.6| ) with s = —2 and s = +2 respectively which are specified by the 
boundary conditions 
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^m^2g-ia;r,^ as r 2M 



-2RZir) - { , (3.7a) 



and 



AI'P A-^e"^'^^* + e+^'^^- , as r 2M 
+2i?rJ(r) ~ <! . (3.7b) 

5;'Jr-V'^^*, asr^oo 

In the light of ( p.9|) , ( p.l2|) it is clear that _2-R;°(r) and +2-RJ2'('^) uniquely specified by 
the above conditions; explicit formulae for and Bj^ will be produced in a later section. 

Use of the subscripts "in" and "up" derives from the analogy with classical scattering 
demonstrated in the previous section; from (|3.7a|) we can interpret _2-RJJi(r)e~*'^* as a unit- 
amplitude spherical wave propagating inwards from infinity and being partially refiected back 
out to infinity and partially transmitted across the horizon, whereas from (|3.7b|) we see that 
+2-RjJJ;'(r)e~*'^* represents a unit-amplitude spherical wave propagating upwards from the past 
horizon, and being partially refiected back and partially transmitted out to infinity. 
is the refiection coefficient and and |-B/|^P is the transmission coefficient for the scattering 
process. The situation is depicted in Fig. |l|. 

The fact that the "in" and "up" modes (|3.5|) are expressed in different gauges will not 
cause any difficulty later as we shall only use the modes to construct objects which are gauge 
independent (namely |\i/^|^). 

We now verify explicitly that Chrzanowski's mode set comprises a complete set of so- 
lutions to the perturbed field equations for Schwarzschild space-time, and in the process 
derive a set of equations which will prove valuable for calculating expectation values of the 
perturbed Weyl scalars later. 

First substitute ( |3.5aD for /i(a)(b) in Eq. ( |3.3aD . Since we are working in the ingoing 
radiation gauge for which /i(a)(i) = 0, we have simply 

^0 [h'^l, m, uj, P- x)] = ^{D-p){D- p)/ig)(3) (/, m, oo, P; x) 

= -^N'^iD -p){D- p){D -p){D + 3p) +,Yr{e, 0) -2i?;°(r)e— * 

where ( |3.5a| ) and ( |3.6| ) have been used in obtaining the second line. 
Recalling that 

^ „r^ r'^ d d 

and p = —r~^, we see that Dp = p^ and consequently 

1 



h'^{l,m,uj,P;x) 



--N'^DDDD +2Yr -2i?;° (r)e-^"* 
-iiV- +2Yne, (l))e-'-'VVVV _2i?;° (r) (3.8) 



where 
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V 



dr 



luor 



A 



(3.9) 



We now use the following result of Press and Teukolsky; [jTSl if _2-Ri(^(r) is any solution of 
the radial equation ( p.6| ) with s = —2, then WW -2Riu}{j) will be a solution of ( |2.6| ) with 
s = +2. So let +2-R/Ji('") t>e that particular solution of the s = +2 radial equation which is 
given by 



AWW^2RZ{r) 



(3.10) 



(the 4 is present to achieve consistency with the normalizations used in Refs. |15] and 
Then (|3.8| ) becomes simply 



\i>o h''^{l,m,uj,P;x) 



iiv^%2i?;:i(r)+2i^r( 



(3.11) 



In a similar manner one can substitute (|3.5a|) into Eq. ( p.3b|) for ^^4, and also (working 
in the outgoing gauge now) insert equation ( |3.5b| ) for the "up" mode into both ( |3.3a] ) and 
( p.3b|) . Altogether one obtains 



^T{l,m,uj,P-x) 



-iujt 



_ ATin (T^ec*,. + 12^M.;P) ^^K^ir) ^2Yl 

o 



-iuit 



-iujt 



where we are using the obvious notation m,u,P; x) for h^{l,m,u!, P; x) 

= (/-l)/(/ + l)(/ + 2) + 122Mcu . 
_2-R"^(r) is defined by the equation 



^R^Sir) 



4|C, 



lu) 



up / 



(3.12a) 
(3.12b) 
(3.12c) 
(3.12d) 

and 
(3.13) 

(3.14) 



and is known [|T^ to be a solution of the s = —2 radial equation. We stress that since \E'^ 



is gauge independent, the final expressions ( |3.12| ) do not depend on the particular gauges 
chosen to facilitate their derivation. 

Eq.( |3.12| ) is clearly in perfect agreement with Teukolsky's solution (^.5]), and we conclude 
that the mode set ( p.5|) does indeed constitute a complete set of solutions to the equation 
of motion for the metric perturbation. 

The constants A^^ appearing in ( |3.5| ) are fixed in the quantum theory by demanding that 
the mode set satisfy the othonormality conditions 

( }— h^{l, m, u, P] x), }— h^'{l', m', u', P'; x)) = 6AA>6ii>6mm'S{uj - uj')6pp> (3.15) 

ylOTT VlOTT 



where the inner product ( , ) is defined as follows for arbitrary complex symmetric tensor 
fields ipai3 and (j)af3- 

^) = ll {r^*V,4>ap - + 24>a,Vprf'* - 2C^V;3r^) (3.16) 

(V'q/S) (pai3 are the trace- free parts of ipaf3, <Pai3 respectively). S is an arbitrary spacehke 
hypersurface in the exterior Schwarzschild space-time, so one can divide the left hand side 
of ( p.l5|) into three integrals 



H- 

where S is the surface enclosing the volume V shown in Fig. |^. 

The integral on the closed hypersurface S converts to an integral over the space-time 
region V via Stokes' theorem; this subsequently vanishes because its integrand is identically 
zero by virtue of the linearized field equations (to which h^j^{l,m,uj, P; x) are solutions). 
The left hand side of ( p.lSj) thus splits conveniently into just two integrals over TC~ and X~ , 



each of which can be evaluated explicitly since one knows the form of the radial functions 
_2-R|Ji, +2-Rz^ at both the horizon and spatial infinity (see (|3.7|) ). To proceed rigorously one 
transforms {t, r) to Kruskal null co-ordinates 

U = _Q{r,-t)/{iM) ^ y ^ g(r.+i)/(4M) 

and employs 

4Mr^ 

(dE'^)^_ = -6(2)^— ^d?7sinMM0 

(dS^')i- = e(i)'^^^dV^sin^dM0 

as future-directed surface elements for the null hypersurfaces X~ and 7i~ respectively. In- 
tegrations over angular co-ordinates can be performed using the orthonormality relations 
( [B4D for the spherical harmonics. One discovers that Chrzanowski's mode set satisfies the 
orthonormality conditions ( |3.15|) provided 



2 



2 16 



(3.17a) 



|iV"Pr=- —— (3.17b) 

where 

= 2Muj (l + AM'^uj'^) (l + IGM^cj^) . (3.18) 

These normalization factors differ from those given in Ref. @] only in that we are working 
in units where G = 1 while the authors of Ref. chose units where IGvrG = 1. 

An arbitrary linear perturbation of the Schwarzschild background can be expanded as 
follows in terms of the complex mode set (|3.5| ): 
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K 



(3.19) 



where is a shorthand for (A, /, m, tu, P). When the theory is quantized h^^, ax and a|f 
become operators on the Hilbert space of quantum states of the system, subject to canonical 
commutation relations which (since the mode set is orthonormal) take the simple form 



(3.20) 



Consider now the expectation value where \B) is the vacuum state associated 

with the mode set ( |3.5| ), called the Boulware vacuum. Since \E'a is linear in /i^i, and its 
derivatives (see (|3.3| )), one can substitute expansion (|3.19|) for /i^,^ and obtain 



{B\\^Af\B) = {B\Y,{aK^A 

K 

K' 

= E [h 



} 



X 

K'* 



h 



}\B) 



Written out fully this is 

{B\\^Af\B) = Y.Y. / {\nil,m,u;,P;x)f + \^J{l,m,uj,P;x)f} (3.21) 

l=2m=-r^ P=±l 

where "^^il, m,Lj,P; x) for A = in, up are given explicitly by ( p.l2| ). 

Although the Boulware state \B) corresponds to the usual notion of a vacuum at large 
radii, it becomes pathological as one approaches the horizon, in the sense that the compo- 
nents of the renormalized expectation value of the energy-momentum tensor for the scalar 
and electromagnetic fields in the Boulware vacuum diverge in a freely-falling frame as 
r ^ 2M (Ref. 



see also 17 , 18 



I). It is thus only appropriate for studying a body 
whose radius exceeds 2M such as a neutron star. 

The Unruh vacuum \U) is the vacuum state of a mode set whose "in" modes are positive 
frequency with respect to the co-ordinate t (i.e. they are precisely the /i™j,(Z, m, u;, P; x) 
defined by ( p.5a]) ), but whose "up" modes are defined to be positive frequency with respect 
to the Kruskal null co-ordinate U. The renormalized energy-momentum tensor for the scalar 
field in this state is regular on 7i+ but not on T-C~, and corresponds to a fiux of black body 
radiation as r — oo. The Unruh vacuum approximates the state of the field long after the 
gravitational collapse of a massive body. 

Finally the Hartle-Hawking vacuum \H), defined to be that of a mode set whose "in" 
modes are positive frequency with respect to V and whose "up" modes are positive frequency 
with respect to U, yields a renormalized energy-momentum tensor which is regular on both 
horizons, but one then has a bath of thermal "Hawking" radiation at infinity, so that \H) 
is not a true vacuum state in the usual sense. It corresponds instead to a black hole in 
unstable equilibrium with an infinite bath of thermal radiation. 

If the calculation leading to (|3.21|) is performed using the Unruh and Hartle-Hawking 
mode sets instead of the Boulware mode set (13.51), one obtains [HI 
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oo +1 



{u\\^Am=Y: E / E 



and 



1=2 m=-l 



X 



P=±l 



[\¥X{l,m,uj,P]x)f + coth{4nMLj) m, tu, P; (3.22) 



oo +Z 



1=2 m=-l 



H\\mA\ \H)=Y. J2 coth(47rMcj) 



p=±i 



X 



{ m, ^, P; x) f + l^Jil, m, u, P; x) |'} 



(3.23) 



respectively. On account of their distributional character the product of two field operators 
^(a)(fe) (oi' their derivatives) evaluated at the same space-time point is a priori ill-defined; 
as a result the expectation value of any such product is infinite. Therefore, since \'^a\ 
is quadratic in /i^j, and its derivatives (see (^]3|)), expressions (|3.21| - p^23D are all infinite. 
Rather than become imbrued in the construction of a renormalization scheme for 
analogous to that employed in the evaluation of {T^'^) for quantum fields of lower spin, we 
elect instead to work with the differences 



POO 

2E E / do. E 



1=2 m=-l 



P=±l 



(3.24a) 



and 



2\U-B. 



{U\\^a\'\U)-{B\\^a\'\B) 

oo +1 

2E E / E 



1=2 m=-l 



^Jil,m,LU,P;x) 

(oSnMuj _ 1 ^ 



P=±l 



(3.24b) 



These differences are automatically finite for r > 2M, by virtue of the fact that the divergent 
component of the expectation value of such a two-point product is purely geometrical and 
hence state-independent. 

Using these differences we can, of course, work out the difference 



^^1')^-^ = {H\\^a\'\H) - {B\\^a\'\B) = {I^aH 



2\H~U 



2sU-B 



It is a consequence of the time-reversal invariance of this difference and the transformation 
properties of the Kinnersley tetrad under time-reversal that 



{r-2Mr{\^on 



2\H~B 



16r^(|^4r) 



2\H-B 



(3.25) 



Substituting (|T|) for ,m,LJ, P; x) in (|3.24|) we obtain finally 

2 

(l^olY 



1 oo 



UJ^ {fS^Mu^ - 1) 



(3.26a) 
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(1^410 



2\H-U 



-| oo 



duj \Ci 



(eS^Mc. _ l^ 

2 



1 oo „ 



(3.26b) 
(3.26c) 
(3.26d) 



where the addition theorem ( p6|) for spin-weighted spherical harmonics has been employed 
to perform the sum over m, and ( |3.17| ) has been substituted for |A^^p. 



IV. POWER SERIES REPRESENTATION FOR sRiu 



In order to evaluate the expressions (|3.26|) , an explicit representation of the radial func- 
tion sRiui valid throughout the region exterior to the black hole is required. In what follows 



we use the methods of Leaver |]T9[ to construct a general solution of the spin s radial equa- 
tion (|2.6| ) in terms of a convergent power series about r = 2M. This procedure has much 
better stability properties than performing a straight numerical integration of the differential 
equation (|2.6| ). 

In light of ( |2.12| ) one particular solution of the spin s radial equation (|2.6| ) can be ex- 
panded as 



where 



k=0 



(4.1) 



(4.2) 



The coefficients ak{l, uj, s) are determined by substituting ([4.1|) into the radial equation ( p.6|) 
and changing variable from r to x = 1 — (2M/r). This yields the following differential 
equation for sSiuj'. 



X 



1 - + {(1 - s) (l - x^) - 2xil - x) - AiMu] 



+ i{2s~l)^^ + s{s-l)-l{l + l) 



dx 

sSluj= 0. 



(4.3) 



Substitution of expansion ( |4.1| ) for gSi^ then yields the following three term recursion relation 
which determines the ak{l,uj, s): 



k{k — AiMuj — s)afc 
+ [-3(A; - i f + {k-l){s + AiMuj) + AiMuj{2s - 1) + s{s - 1) - /(/ + 1) 



3(A; - 2f + s{k - 2) - s{s - 1) + /(/ + 1) 



a'k-2 + 



_(A;-3)2-s(A;-3) a^-s = (4.4) 



with the initial conditions = 1, a_i = a_2 = 0. 



13 



Another independent solution of the Teukolsky equation is (c.f. Eq. ( p.l2D ) 

MUr) = e+'"'-* sSUr) (4.5) 

with 

sSUr) = Y.~a,,{l,uj,s)[l (4.6) 

fc=0 ^ 

say. Substitution of ( [4.5| ) into ( p.6| ) yields the complex conjugate of ( [4.3| ) with s replaced by 
— s, so that 

sSiuj = -sSi^ . (4.7) 
One concludes that the general solution of the spin s radial equation ( p.6|) is 

MUr) = ci A-^e-'-^* ,SUr) + C2e+^-'-* _,S;jr) (4.8) 



where sSiuj{r) is given by ( ^TTl) , the coefficients ak{l,uj, s) being determined by the recursion 
relation (|]|I). 

The choice of coefficients Ci, C2 which yields either _2-RJ^(r) or +2-R"^(r) is determined 
by comparing the form of either of these particular solutions as r — > 2M (see Eq. ( |3.7| )) with 
that of the general solution (WM). One deduces the representations 



RZ= 5- A^e— _2^^. (4.9a) 



-2^1, 

+2RZ= AL'A-^e— +2^^. + e+^-- -A (4.9b) 

which are valid throughout the exterior region. 

The corresponding representations of +2-RS and _2-RJaj' are obtained by substituting 
( 4.9a ) and ( 4.9b| ) into ( 3.10| ) and ( p.l4| ) respectively, and performing the derivative operations 
explicitly; derivatives of gSi^^ of order > 2 can be eliminated with the aid of ( [4.3[ ). The result 
is 

+2i?;°= 4S|:iA-2e-^-'^* {ai^ ^^Si^ + A.A ^^Sl} (4.9c) 

+e+^-- {7/.A _25L' + -2SI} (4.9d) 
where the prime denotes d/dr and 

ai^{r) = (/ - !)/(/ + !)(/ + 2)A2 - 2iujr'^ {(/ -!)(/ + 2)r (At^ - lOMr + 47^^) 

+10Mr2 - 40MV + 24M3} - 4cuV {lOMr - 24M2 + 3(/ - !)(/ + 2)A} 
+Siuj^r^{r + 2M} + 16cjV^ 

/?;^(r) = -2icjr2 (4Mr - 12M2 + 2{l -!)(/ + 2)A) + ^iuh^ 

7;^(r) = 2iujr'^ {4 (2r'^ - 3Mr - 3M^) + 2{l - 2)(/ + 3)A} - 8iu^r^ 

Si^(r) = {l- !)/(/ + !)(/ + 2)A - Aicur {Sr^ - 9Mr - GM^ + (/ - 2)(/ + 3)r(2r - 3M)} 
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-Auj^r^ {2(2r + 3M) + (/ - 2)(/ + 3)r} + 2Wr^ 

+2iur^ {2 (Sr^ - ISMr^ - 12MV + 12M^) + 2(/ - 2)(/ + 3)(2r - M)A} 
-4cjV^ {2 (er^ - 7Mr - 12M^) + 3(/ - 2)(/ + 3)A} 
-Siu^r^r + 2M) + 16tuV^ . 

It is reassuring to check that in the hmit as r — > 2M these latter representations reduce 



to 



pm 



1«(2M)V , 

{l + 2iMu) 
i{l + 2iMuj) 



16(2M)4p^ 



\Ci^\' {l-2iMu) 



(4.10a) 
(4.10b) 



which is consistent with their expected behavior (|2.12| ) in this hmit. Their behavior as 
r — ^ cxD is derived in Appendix 0. For completeness we give the leading behavior here 



+ 



p 12 p+iojr* 



(4.10c) 
(4.10d) 



One can now insert the above power series representations ( ^.9| ) for the radial functions 
into the expressions (p.26|) for the differences in expectation values of the perturbed Weyl 
scalars, thereby obtaining 



1=2 
/\ 4 00 



^5 CgSTrMo; _ ^ 



/■oo rl/ 1 I Ri'^ |2 



2^nr^ 

i2^u-B^ 1 ^^ y(2/ + i) /' 



1*410 



2\U-B 



\Al^\'A-^\+2SiJ'+\-2SiJ' +27^e 



(gSTrAfo; _ I 

lui 



^;„.e LA +2^luj ~2^liu 



(4.11a) 
(4.11b) 

} (4.11c) 



1 ^ roo du; 

2Hr^{2My f^} ^ 'Jo (eS-^^- - 1) 



+27^e 



■^L"e-2^"^*A-M7z.+25L + 5/.+25.J(7LA-25L + eL-25,j]} . (4.11d) 



V. THE REFLECTION AND TRANSMISSION COEFFICIENTS 



It remains to derive formulae from which Bf^ and can be determined explicitly; 
expressions ( [4.11| ) for the expectation values can then be evaluated numerically. 
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Dm 



Such a formula for Bf^ is derived by comparing the power series representation for _2 
(which is vahd for all r > 2M) with its asymptotic expansion valid only at large radii 
(asymptotic expansions of the various radial functions are derived in Appendix 0). Thus 
from (|4.9a|) and ( |Cl|) we have at large r 



+ 



1 ei 62 63 64 
- + — + — + — + — + 

tyt ty2 iy^3 jy^ 5 \ lyy. 



(5.1) 



with 



ei= bi + 2(2M) 

62= 62 + 2(2M)6i + 3(2M)2 

63= 63 + 2(2M)62 + 3(2M)26i + 4(2M)^ 

64= h + 2(2M)63 + 3(2M)262 + 4(2M)^6i + 5(2M)^ 

where the bi are given in Appendix 0. By ignoring terms of order on the right hand 
side of this equation one could write down an approximate formula for the ratio A^/B]^, 
evaluate it numerically at some suitably large value of r, and hence determine |-B;™P from 
the Wronskian relation 



75 m 



26cj3 



Am 



(5.2) 



this last equation being a variant on one of a complete set of Wronskian relations between 
the reflection and transmission coefficients which are derived in Appendix 0. However a 
more accurate approximation can be obtained by first applying the operator (recall 
(|]9D) to both sides of (|]l|), so that one need only ignore terms of order r ^ to arrive at the 
approximation 



Am 

B. 



liU 



/l_|_^_|_^_|_/4_|_^ 



-2iu)r, 



dr 



(5.3) 



for the ratio of the incoming coefficients, where 
/i= 22t^ 

/2= 2iuo [61 + (2M)] - 1 

^62 + (2M)6i + {2Mf] - 261 



/4= 2iuj 



63 + (2M)62 + (2M)26i + {2Mf] - 362 

64 + (2M)63 + (2M)262 + (2M)=^6i + {2MY 



46. 



In practice the right hand side of ( |5.3| ) is evaluated for large and increasing values of r until 
it has converged to the desired accuracy; the result is then inserted in (|5.2| ) to yield |-B;°p. 
A slightly different procedure is used to evaluate the reflection amplitude A^^. A com- 



parison of the power series expansion (|4.9b| ) for +2-R/^ with its asymptotic form 
large radii yields 



at 



16 



+20iui + A e -20i^ - e 



tj-t .y^ 2 ^4 .y^5 ^ 



with 



g2= ci - AM 
g^=C2- 4Mci + 4M^ 
^4= C3 - 4MC2 + 4M2ci 
^5= C4 - 4MC3 + 4M2c2 

where the are given in Appendix y, which when operated on by T> gives 



( +25'ic 

dr 



2iujr 
A 



+ Ae^ 



A-(-25rJ+4(r-M)_25r^ 



To this last equation one now adds 



/^O lyyj 

multiphed by 



ui \ ^ h2 h-i hi 

ijn ryj ryj 



where 



(5.4) 



(5.5) 



h2= 92 

h4= 35(4 - g3h2 - g2h 

h5= 45-5 - g4h2 - gsh - g2h4 

so that the right hand side of the resulting equation is O (r~^) and can be ignored at large 
radii. An accurate approximation for is therefore given by 



/I up 



-A'e 



2'2iuT* 



dr 



( -2SI) + h{r) + 



X 



^(+2^,.)+ Mr) 



4(r - M) 
A 



-2S, 



2^luj 

1 



A 



(5.6) 



which is evaluated in the same way as (|5.3| ). 

The values of the reflection and transmission amplitudes, which we have computed nu- 
merically, are displayed graphically in Figs. |^ through ^. We have verified that these values 
are correct using the following (independent) checks: 

(i) The rate of decrease of mass of the black hole due to graviton emission is given by the 
formula lOll 



M 



(5.7) 



where 
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1 



dE^°'/dt 



oo +1 



(5.8) 



Substituting first (p.l5|) , (p. 16a]) for tlie energy fluxes, and then \i/J['(/, m, cu, P; x) (see Eq. 
(CT)) for ^ 0, one obtains tlie following expression for the fractional absorption of incoming 
radiation by the black hole (spherical harmonics are eliminated with the aid of (p^ )): 



dE^°7dt 


1 4u;2 




r=2M 


dPi°/dt 


(2M)4 (16M2u;2 + 




r— >oo 



(2M)' 



where the last line follows from ( [4.10a[ ) and ( |(J2[ ), and hence 

oo 

M., = ^^^^(2/+ 1) 



(2M)5^^^, , pjEiSP 



47r 



«=2 



(e 



1) 



(5.9) 



Using the values obtained numerically for |-B;™| in the above equations one obtains the 
luminosity spectrum shown in Fig. |^ and the value M = (3-84 x 10~^)M~'^ for the total 
luminosity due to graviton emission; both are in close agreement with the previous results of 
Page IP] (Page has M = (3-81 x 10~^)M~2; we believe the slight discrepancy to be caused by 
his choice of an upper limit for the u integral which is not high enough to ensure 3 significant 
figure accuracy.) (ii) When 2Muj ^ 1 the spin s radial equation can be solved analytically in 
two overlapping regions, one of which includes the horizon, and the other of which extends to 
infinity; solutions are expressed in terms of ordinary and confluent hypergeometric functions 
in the respective regions (see Ref. |]2D|). By matching these general solutions in the region 
of overlap, and then specialising to the particular solutions -2-R;Ji, +2Rf^ in turn, one can 
derive analytic approximations for the reflection and transmission amplitudes which are 
valid in the limit as 2Muj —>■ 0. The derivation is performed in detail in Appendix E. The 
exact numerical values for the amplitudes are consistent with the analytic approximation in 
the limit 2Muj — >■ 0. (iii) A final check on the numerical values we have obtained for 
and I^JI^I^ is provided by the Wronskian relations derived in Appendix D. From ( P4|) and 
the following relation between the 'in' and 'up' reflection coefficients may be obtained: 



/I in 1 2 _ IQuj''^ I . up 1 2 



(2M)9(16M2cj2 + l)p, 

Substituting this relation into (P5|) then yields 



(5.10) 



16(2M)8p2 I zo^l + 



I Din 1 2 



\Bl 



(5.11) 



We have shown our numerical results to be consistent with the above equation. 
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VI. ASYMPTOTIC ANALYSIS 



Before we discuss the numerical evaluation of Eqs. (|3.26|) , it is useful to consider their 
asymptotic behavior in the limits r 2M and r — » oo. 

Consider first the limit r — > 2M. The leading behavior of the various radial functions in 
the limit as r — 2M may be deduced from Eqs ( [4 .91 ), and ( [4.10| ); Eqs. ( |3.26| ) then give 



9 roo (^ _L 1 R A/r2,.,2\ ^ lRm|2 



(|$„|^)»-<'.li^i^V(2;+l) rd.'^ + "'P"lf-l (6.1a) 

^' ^ ttA^ 'Jo cu^ (eS^A^"^ - 1) ^ ' 

A 4 oo 1/^ |2|r5in|2 



The integrals in Eq. ( |6.1|) may be evaluated numerically to give: 

|2\/f-l/ 



lirn^ (1 - 2M/r)^(|^or)-^"^~ 1-96 x 10"^ (2M)-6 (6.2a) 



lirn^ (1 - r/2M)-4(|^4r)^-^^ 9-14 x 10"^ {2My (6.2b) 

The same method cannot be used to evaluate the limit of as r 2M. For 
writing 

oo 

(I^a|') = E(2^ + 1)^az (6.3) 

1=2 



one finds (by substituting (|C^ ) in (|3.26a)) ) that 

.1(271/ 
2^157? 

and 



In neither case will the sum over / converge since the terms approach a constant times 
(2/ + 1) for large I (note from (|6.4a| ) that we have already seen that the corresponding sum 
over ipQi~^ converges). The problem is that the limit as r — >■ 2M and the sum over I do 
not commute in this case since our asymptotic expansions for the Teukolsky functions are 
not uniform in /. Nevertheless Eqs. ( p. 41) do provide a useful check against the numerical 
results. 

To find the asymptotic behavior of (|\i/o|^)^~'^ and (|\i/4|^)^~^ as r ^ 2M we need 
asymptotic expressions for the Teukolsky functions which are uniform in I. To this end we 
follow Candelas [0. Near r = 2M the Teukolsky equation (|2.6|) may be approximated by 
the equation 



d^ ^ d 2Muji2Muo -is) il-sMl + s + l) 

,_2M)- +(.+!)-+ --^i — - 

Defining ^ = (r/2M — 1)^^^ this becomes 



sRUr) = 0. (6.5) 
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+ 



{2s + 1) d 8Muj{2Muj 



+ 



IS) 



-4(/-s)(/ + s + l) 



sRiJr) 



(6.6) 



which admits solutions in terms of modified Bessel functions. Since it is clear from Eq. 



that the asymptotic forms of (|^E'o| ) and (|^'4| )^~^ as r — 2M are determined by the 
contribution from large /, we may further approximate the constant term in the potential in 
Eq. (|6.6|) by —AP and then the solutions are given by 



s+4iMuj 



and 



These solutions are uniformly valid for large /. 
We now concentrate on (|\i/o|^)^~'^ and write 



~ air^K2+4iMu.m) + f^ir^i-^-^iMMio- 



Following the arguments of Candelas, as / ^ oo for fixed C,, +2-Rj^(0 ~^ and so Pi is an 
exponentially small function of / that will not therefore contribute to the leading asymptotic 



behavior of (l^&ol ) 
as e — > yields 



as r 



ai 



— > 2M. Furthermore comparison of (|6.7| ) with ( p.7b| ) for fixed I 
2iT{3 + AiMu) sinh(47rMcj)e2**^'^ 



(6.8) 



The elementary identity 



then enables us to write 



Since for large /, 

oo 



|r(3 + 4iMcu)|^ 



5m 



\ai 



sinh(47rMu;) 
2 _ 32p<^sinh(47rMcj) 



(6.9) 



it follows that to leading order 
64p^^ sinh(47rMcj) 



1=2 



(6.10) 



where we have used Eq. (6.576.3) of [^. Inserting this into Eq. ( ^.26cD we obtain the 
leading asymptotic form as r ^ 2M as 



A'(i^or) 



2\U-B 



16(2M)' 



duj- 



57f{r - 2M) Jo [e^^Mu; 

191 (2M)3 
24 3157r (r - 2M) 

2 mf 



(6.11) 



1-21 X 10" 



(r - 2M) 
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or equivalently, 



lim (r/2M- l)^(l^or)^"^ ~ 1-21 X 10-^ (2M)"^ (6.12) 

r^2M 



Since, by (B.la| ), A^(|\&o|^)'^ ^ is finite on the horizon it follows that the asymptotic form 
of A^d^ol r^^^ as r ^ 2M is also given by (|6lll ). Then from 



,2,H-B 191 1 



I^J^)^-^ ~ (6.13) 

' ^' ' 28 3157r (2M)5(r-2M) ^ ' 



7-54 X 10"^ 



1 



(2M)5(r-2M)' 



Since, by ( |6.1b[ ), (|\i'4|^)^ ^ vanishes at the horizon it follows that the asymptotic form of 
(|^^|2)C/-B as r ^ 2M is also given by (|67[3| ). 

We now turn to the asymptotic forms at infinity. Substituting the asymptotic expansions 
(^D and (|C|) for +2RZ and .2RZ respectively in Eqs. (jOe^ ), ( [OMD yields the following 
formulae which describe the leading behavior of {\'^ a\^)^~^ in the limit as r — 00: 

1 00 1/^ |2| D^iPlS 

(I'^oI'>"'°~ „-,oL... I:P' + 1)/ d., !^'" (6.14a) 



47r(2M)Vio ^ Jo {q&^Mu> _ 

Using the Wronskian relation (p4|) , the last equation may be rewritten as 



(1*^^"-" ~ ^ £(2' + 1) r d^ ..,g:iS'\, (6.14c) 



Eq. ( |6.14c| ) is in accord with Eq. ( |5.9| ) for the luminosity since Eq. ( ^.16b| ) implies that for 

r ^ 00 

(1^4.1^)^-^ - (6.15) 
The integrals in Eq. ( |6.14|) may be evaluated numerically to give: 

Jiim {r/2MyW^o\Y'^- 2-91 x 10^^ {2M)~^ (6.16a) 
Jiim (r/2M)2(|^4|2)^-^^ 8-42 x 10"^ {2My^ (6.16b) 

Finally, we consider the limit as r — »• 00 of 

ii^MY-'' = ii^MY-'' + {\^A\y-^. 

We find that that 

^^^-^ ~ 16^1^-^ ~ ^, (6.17) 

•^0/ ^41 24(2M)n57rr2' ^ ^ 
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where Wronskian relations ( |D4|) and (P5| ) have been used. As before these forms provide 



a useful check on our numerical results but are not good enough to yield the asymptotic 
forms. 

The simplest way to obtain the asymptotic forms at infinity is to use the fact that at 
infinity the calculation reduces to a problem in Minkowski space. Then quantum effects are 
negligible (except in providing the appropriate temperature) and we may tak 

e over the standard classical results given by Eq. ( |2.16| ). Firstly, using spherical symme- 
try, we have 

l^o..!' = -647iu;^T^\' and |^4.|' = 47iu^T^^\' 

where T™^* and T°"*r* denote the energy flux components of frequency u in the energy- 
momentum tensors associated with the incoming and outgoing gravitational waves, respec- 
tively. Then, since we are dealing with a thermal bath of gravitons at the Hawking temper- 
ature, l/(87rM), we have 



-(7;^V) = (T-V)=2x J 



2 X 



1 UJ^ 



where the factor of 2 arises from the number of polarization states, and multiplies the ingoing 
and outgoing energy flux for a real massless scalar field in a thermal state with the Hawking 
temperature in Minkowski space. Thus, as r ^ oo we have 

(1^0^)--- 16(1^4^)-- 1^ fd. ^^j'_^^ (6.18a) 
^ {2M)-^ (6.18b) 



20167r 

1-58 X 10"^ {2M)-^ (6.18c) 



Since, by Eq. ( |67HD (|^o|^)^~-^ and (|^4|^)^~-^ vanish at infinity, Eq. (^JEj ) also gives the 
asymptotic form for (|\i/o|^)^~'^ and {{"^if)^'^. 



VII. NUMERICAL RESULTS 



In this section we present the results of evaluating expressions (^4.11|) for the differ 



2 

ences between renormalized expectation values of \'^a\ in the Boulware, Unruh and Hartle- 
Hawking vacua, using values for the coefficients 1-8;™^ and A^^ which were obtained numer- 
ically in Sec. 0. 

All power series and sums over / converge swiftly; likewise the infinite upper limit on 
the integral presents no difficulty, and may be approximated with sufficient accuracy by 
(2M)a; = 2 in every case. There are nevertheless technical problems associated with the 
computation, two of which merit some explanation here. 

The first problem concerns the evaluation of the integrands of ( [4.11| ) at the lower limit 



of the integral uu = 0: although the integrand of ^ can be seen to vanish at = 
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(consider ( |E7|) ) , the integrand of {{"^ a^)^'^ is non-zero at = and its value there can 
only be determined by expanding the integrand as a Taylor series in powers of u [0. To 
this end consider the Taylor expansions 

-2Siu;= -2S10 + [~2Sm]'uj + O (u;^) 



4UP 1 



LU 



where we have introduced the notation 

[fio]' - 



dfi 



duj 



(7.1) 



J a;=0 



We have expanded +2S1UJ rather than +25';^ since the latter diverges like uj~^ as u —>■ 
(consider Eq. (|4.4|) when k = s = 2). We have also used A^q = which follows from (^). 
Inserting these expansions into the power series representation ( |4.9b| ) for +2Rf^ one obtains 



UJ 



In order that Eq. (|3.26c|) for (|\i/o| )^ be finite for finite r > 2M, we must have +2R1UJ 
as u where n > 1, i.e. 



(7.2) 



^2R7 = 



and the value of the integrand in expression ( p.26c|) for (|\l/o| ) at = is 

1 _,_.2 

2%^(2M)' 



-^E(2/ + l)[(/-l)/(/ + l)(/ + 2)f [+2i?jr 



(7.3) 



(7.4) 



where [+2-R"o'']' is the coefficient multiplying u in (|7.2|). ( |7.4| ) can now be evaluated numeri- 
cally; in particular the power series are determined straightforwardly from (|4.2|) and ([4.4|), 



and [A^q]' is obtained from the formula 



[A 



upy 

10 J 



- -2'5';o 
+2>S';o 



(7.5) 



which follows from (^) and Q. The second derivative [Af^]" = [{dyd^cu) {AZ)]^^o is 
computed using the approximation 



[A 



10 . 



+2S10 



d 



d 



+ ( -2SU + 2tn -2SL 



(7.6) 



1^=0 



(which follows from expanding both sides of ( |5.6| ) in powers of u and using ( [7. 5| ) ) . Note that 
( [r.6|) will only become independent of r at large radii, in contrast with ( [7.5|) which may be 
computed at any value of r > 2M. 
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The value of the integrand of ( p.26d| ) at = is found in a similar way, requiring 
only that (d/dr)_|_2'S';a; and [(d/dr)_25'f^]^=o evaluated numerically in addition to 
the power series considered above. In both cases the result is finite and non-vanishing. 

The second technical difficulty also arises during the computation of when 
one attempts to evaluate | s-R/^P numerically using the power series representations ( [4.9bD , 
( [4.9d|) . This latter quantity is composed of two pieces (see e.g. the expression for | +2-R"^P 
in braces in (|4.9c|) ), which turn out to be opposite in sign but equal in magnitude to high 
accuracy; one must therefore work to high precision in order to produce reliable results. 
The problem becomes more severe as — > 0. This difficulty can be understood by once 
again invoking the analogy between our system and a classical scattering problem; at low 
frequencies, upcoming radiation from the black hole is unable to surmount the potential 
barrier sV/t^(r) and is instead completely reflected back across the event horizon. 

The flnal results are displayed graphically in Figs. ^ |ll], where the quantities have been 
scaled to give flnite values on the horizon and at inflnity. The most physically interest- 
ing graphs are those for {\'^4f)^~^ (Fig- 0) corresponding to outgoing radiation from an 
evaporating black hole and Af)^~^ (Fig- 0) corresponding to a black hole in thermal 
equilibrium at the Hawking temperature. The most striking feature of these graphs is the 
very rapid decline in vacuum activity with r; this is even more pronounced when the asymp- 
totic scaling (which softens the effect) is removed. This suggests that quantum gravitational 
effects may play a highly signiflcant role in determining the back-reaction near the horizon 
even though the asymptotic flux to inflnity measured by M is less that that due to lower 
spin flelds. 



VIII. CONCLUSION 

In Sec. we chose to apply Leaver's method directly to the Teukolsky equation, this 
has the great virtue of directness. We now briefly mention an alternative procedure which 
we considered employing but rejected. This alternative is t 

o work with the Regge- Wheeler equation 



dr* 



0, 



(8.1) 



where 



luj 



r -2M (l{l + 1) 6M' 



\.2) 



Eq. (|8.1|) has independent solution and F^^ with asymptotic behavior 



as r 



-iijjrt, _|_ „in„+it<jrt 



2M, 
as r — >• cxD , 



(8.3) 



and 
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lu) 



g+ic.r. ^ aZe-''^''* , as r 2M, 

(8.4) 

^ic^e"'"*'^''* , as r — > cxD. 



Defining the complex differential operator 



'd = 2r(r - 3M + iur"^) + itu^ + r^Ui^, 



.5) 



we can then write the Teukolsky functions in terms of and F^^ as 

= -^'^FZ (8.6a) 

-^^^ = 1^(1 + 2^^^)(1 + 4^Mu;)^i^J (8.6b) 
+2i?;:i = 16u;^A-2^*i^- (8.6c) 
+2i?zl' = ^^(1 + 2^Mc.)(l + 42MC.) A-2^*F;J (8.6d) 

The apparent advantage of the Regge- Wheeler approach is that it deals with a real 
equation similar to the spin-0 equation. However, the three term recursion relation obtained 
through Leaver's approach is still complex so no real advantage is accrued. In addition, the 
expressions for the Weyl scalars involve derivatives of the Regge- Wheeler functions which 
considerably complicates and obscures the asymptotic analysis required as compared with 
that employed in the more direct approach we have followed. 

Having argued that there is no advantage in using the Regge- Wheeler formalism for our 
calculation, we should point out that some of the asymptotic formulae involving Teukolsky 
transmission and reflection coefficients look far neater in terms of Regge- Wheeler transmis- 
sion and reflection coefficients. For example, it can be shown that 

,2 (2M)5p 4u;5 



and so, Eq. (|5.9|) takes the form 

1 00 / ! I/) P 

M„^2x-g(2,.l,^^Jj^ , (8.8) 

which parallels the scalar result except that the sum starts at / = 2, the lowest radiative 
mode of the graviton field, and there is the extra factor of 2 arising from the number of 
polarisation states for the graviton. 

One might view the research presented here as a precursor to the numerical evaluation 
of the renormalized effective energy-momentum tensor for quantized linear gravitational 
perturbations of a black hole. However, there are severe problems of gauge invariance in 
defining such an object. One approach to this problem, using the Vilkovisky-DeWitt off- 



shell gauge- invariant effective action, was suggested in Ref. [|22[. Since both the Vilkovisky- 
DeWitt effective energy-momentum tensor {T^'^) and {{"^aI ) consist of terms which are 
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quadratic in the metric perturbation h^i, and its derivatives (compare Eqs. ( |3.3|) of this 
paper with Eq. (2.6) of Ref. [^|), the evaluation of the differences in the renormahzed 



expectation values {T^^)^~^ and {T^")^~^ will be a straightforward (though laborious) 
extension of the calculation outlined in this paper. However, if one wishes to compute 
an individual renormalized expectation value {{H\T^'^\H) n say) then new ground must be 
broken, since in this case explicit renormalization is necessary. Problems then arise since the 



graviton renormalization scheme of Ref. can be only implemented in deDonder gauge 
h^v■,^l — \h^^l■,u = (where the propagator has Hadamard form), whereas a complete set of 
solutions to the linearized field equations currently exists only in radiation gauge (see Eqs. 
(p.5|)). If these technical problems can be resolved, it would be interesting to compare this 
quantity with its scalar and electromagnetic analogues. To date it has only been possible to 
make quantitative comparisions of the black hole luminosity of gravitons with that due to 
radiation of massless particles of lower spin |^ . One finds that the luminosity for gravitons 
is less than for lower spin fields but one expects that this is due to the increase in the height 
of the effective potential barrier with spin and so cannot be used to draw conclusions about 
their importance near the black hole. Indeed, we expect that as in the case of quantum 
fields propagating near a conical singularity, the contribution of gravitons will be seen to 
dominate the back-reaction near a spherical singularity (see also the remarks at the end of 

Sec. [yni) . 

In this paper we have chosen to sidestep these gauge-invariance problems and instead, 
following Ref. |^, have concentrated on the gauge- independent Newman- Penrose scalars 
\l/o and \l/4. These scalars are in some sense the true gravitational field variables in the 
Schwarzschild background, for example, it is in principle possible to reformulate the quantum 
theory at the level of the one-loop effective action in terms of '^q and \E'4 (strictly speaking, 
\I/o may be regarded as representating the two radiative degrees of freedom and as acting 
as a 'superpotential' from which \&4 may be obtained). Most importantly, the expectation 
values we have computed provide important physical measures of the vacuum activity of 
the quantized gravitational field around a black hole and, in the asymptotic regimes, they 
directly measure the one-loop quantum gravitational energy flux across the horizon of the 
black hole and to infinity. 



APPENDIX A: NOTATION AND NP CONVENTIONS 

For convenience we repeat here the definition of various symbols used in the text: 

A = r(r - 2M) 
Ci^ = {l- 1)1(1 + i)(/ + 2) + l2iMuj 
= 2Mtu(l + 4:M^iu^){l + WM^co^). 

The following are the conventions and notation employed in this paper for the NP de- 
scription of general relativity. In the NP formalism the geometry of a general space-time 
with metric g^iv{x) is encoded into a null complex tetrad |e(a)^(a;) : a = 1,2,3,4| which 
satisfies the orthonormality conditions 

gij.ue{a)^e{b)' = ri(a){h) ■ 
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V{a){b) is the constant symmetric matrix 



V{a)ib) 




-1 





-10 



1 

1 



(a)(fe) 



and acts as the NP analogue of the metric tensor; in particular tetrad indices are lowered and 
raised by 'r]{a)(b) and its matrix inverse, denoted t]^"-^^''\ respectively. We adopt the convention 
whereby tensor indices are labelled by Greek letters and tetrad indices by Roman letters 
enclosed in parentheses. 



By expressing tensors T..^... in terms of their tetrad components T j^a). 
introducing the spin coefficients (NP analogues of Christoffel symbols) 

7(a)(6)(c) 



and 



e(a)^;7e(6)/3e(c)'^ 



the fundamental equations of general relativity (e.g. Einstein's equation, the Bianchi iden- 
tities) can all be expressed as tetrad equations involving only scalar quantities. 
We also adopt the notation (see e.g. Ref. ||2v 



P = 7(3)(1)(4) 
P = 7(2)(4)(3) 

r = 7(3)(i)(2) 

= 7(2)(4)(1) 



e = ^ (7(2)(1)(1) + 7(3)(4)(1)) 



7=2 W(2)(l)(2) + 7(3)(4)(2) 



e^,fd^ = {D,A,6,6*) 

for the tetrad operators, 

= 7(3)(i)(i) 

= 7(3)(1)(3) 
A = 7(2)(4)(4) 
= 7(2)(4)(2) 

for the spin coefficients, and 

^0 = -C(l)(3)(l)(3) 
^3 = -<^(1)(2)(4)(2) 

for the five independent tetrad components of the Weyl tensor, called the Weyl scalars. 



_ 1 / 

a = 2 (,7(2)(1)(4) + 7(3)(4)(4) 
/3 = - (7(2)(1)(3) + 7(3)(4)(3)) 



^1 = -C(l)(2)(l)(3) 
^4 = -<^(2)(4)(2)(4) 



-C(i 



(1)(3)(4)(2) 



APPENDIX B: SPIN- WEIGHTED SPHERICAL HARMONICS 



Spin-weighted spherical harmonics gYf' 
by the set of equations (Refs fill, [||]) 



are defined for s = —L — I + 1, ... ,1 — 1,1 



s+,Yr{e, 0)= [(/ -s)(i + s + 1)]-- & sYr{0 

s-iYr{9, 4>)=-[{l + s){l-s + 1)]--^ ^' ,Yl^ 



(Bla) 
(Bib) 
(Blc) 
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in terms of the ordinary spherical harmonics 

Here (9" and d'' are operators which act as follows on a quantity 77 of spin-weight s: 



(9" ri= — sin** 9 
d''r]= — sin"* 9 



d id 
-7^ + - 



89 sin 9 d( 
d id 



89 sin 9 ( 



sin 9j 7] 

{{sm'9)r]) 



(B2a) 
(B2b) 



The application oi & on a quantity lowers the spin- weight by 1: the application of on a 
quantity raises the spin- weight by 1. It follows immediately from (|B1D that 



sYr{9, 0) = -{l-s){l + s + 1) sYl' 
or, substituting (p^ ) for & and d'', that 



(B3) 



1 d f . 8' 
sm9— 



1 (9^ 2zscos6' (9 



sin 9 89 y" " 89 J ' sin^ 9 8'^(j) sin^ 61 (90 

-s^ cot^ e + s + (/ - + s + 1) 



0. 



The spin-weighted spherical harmonics satisfy the orthonormality and completeness relations 



(Refs. 0, 0): 



2n 

^0 



d(f) sm9d9sYi"'{9,(j))sYi7'*{9,< 



and 



(B4) 



(B5) 



respectively; these follow essentially by induction on (pl| ) and the corresponding relations for 
the ordinary spherical harmonics. Using the same method one can verify that the "addition 
theorem" 



00 I 

Y: E sYr{9,<i^)sYr{0',<f^')=Si<f)-<i^')5{cos9-cos9') 

l=s m=—l 



2 21 + 1 



m=—l 



4:71 



(B6) 



holds for s = ±1, ±2 as well as s = 0. 25 



APPENDIX C: ASYMPTOTIC EXPANSIONS OF THE RADIAL FUNCTIONS 

AT INFINITY 



Based on ( p.7a|) , will have an asymptotic expansion of the form 
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2RZ - e- 



+ 



+ 6ir^ + fesr + 63 + 



(CI) 



The coefficients tu) are determined by substituting this expansion into the s = — 2 radial 
equation (|2.6| ) and equating the coefficient of each power of r to zero; we find 



(/-l)(/ + 2)z 
2^ 



48a;3 384u;4 
Next by inserting ( PI]) into ( p.lO| ) we deduce that 



Din 



+ 



2 g+jojr* 



-A 



1-^ ^5 



(C2) 



in the hmit as r ^ 00. A similar argument yields the following asymptotic expansions for 
the s = +2 radial functions: 



+2-^-/0; 



(C3) 

(C4) 



as r — s> 00, where 



c,= — [{l-2)il + 3)] 

C2= [{(/ - 3)(/ + 4) - AMuji} ci + 2M {15 - (/ - 2)(/ + 3)}] 



^ 600 



C4-- 



{(/ - 4)(/ + 5) - SMiui} C2 + 2M {30 - (/ - 2)(/ + 3)} ci - 15(2M)^ 



{(/ - 5)(/ + 6) - UMui} C3 + 2M {49 - (/ - 2)(/ + 3)} ca - 24(2M)2ci 



and 



di 



\Cloj\ 

2 



Note that 



2Cii^ 

are consistent with 



d2 
d^ 



l{l + l)iuj 



12Q 



lu) 



APPENDIX D: RELATIONS BETWEEN THE REFLECTION AND 
TRANSMISSION COEFFICIENTS 

Eq (|2.7| ) gives the spin s radial equation in Liouville normal form. The Wronskian 



29 



w 



of any two particular solutions sQiuij), sQiuiij) of the equation will be constant. In addition, 
observe from (^78|) that -sV;*(r) = sViijjij)', thus if sQiujij) is a solution of ( p.7|) , so is 
-s<5L(^)- Hence 

will also be constant. These facts can be used to derive relations between the amplitudes 
< and Bt 

Consider for example 



W 



where we are using the obvious notation 



(Dl) 



This Wronskian is first evaluated on the horizon, by substituting the power series expansions 
(4.7a,d) for -2-^1" and -2-^^" respectively, computing the derivatives, and then taking the 
r 2M limit: we find 



lim W 

r^2M 



:d2) 



where is given by ( p.l8|) . Next we evaluate the same Wronskian in the limit as r — >■ 
oo, using the asymptotic expansions (Pl|), ([C^ ) derived in Appendix |C| for -2-R£, -2-^/^ 
respectively. The result is 



lim W 

r— >oo 



luj 



(D3) 



Constancy of the Wronskian therefore yields the following relation between the 'up' and 'in' 
transmission amplitudes: 



^ pin 



(D4) 



and 



A similar analysis of the Wronskians W _2(5JJi(r), +2QJJi*(r) 
_2Q5Ji,(r), +2Q|JJ**(r) serves to complete the set of independent relations between the 



W 

reflection and transmission coefficients, yielding 



i2uj) 



A 



lu) 



2 ^ (2M) V 



B 



lu) 



and 



'-2zu)Ai:^BZ + 4M (1 - 2iMu^) BZAZ* = 



(D5) 



(D6) 
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respectively. All other relations between reflection and transmission coefficients follow from 
the above three. 

Eq. ( |D5| ) may also be derived from conservation of energy ffux, 

d^i" dE°"* dE^°' 

"dT dt ~ dt ■ ^ ' 



To see this, first substitute ( |3.12|) for m, u, P; x) in (|2.15|) , ( |2.16| ) and integrate over the 



solid angle; next replace ±2-RlJi by formulae describing their leading behavior in the limits 
r — i> 2M, r ^ oo as appropriate; finally substitute into the conservation equation above to 
obtain fpEI). 



APPENDIX E: SMALL lo APPROXIMATIONS FOR ^"p, B}'^ 

In this Appendix expressions are derived for the leading behavior of A^^ and 5;" as u 
tends to zero. Our method is analogous to that employed by Page for his investigation of 
particle emission rates from a Kerr black hole, and in particular makes use of approximate 
hypergeometric solutions of the spin s radial equation obtained by Churilov and Starobin- 
skii . g 



In terms of the dimensionless quantities x = (r/2M — 1) and k = 2Muj , the spin s radial 
equation (|2.(j|) may be approximated when -C 1 by 



d^ d 

x'^{x + + (s + l)x{x + l){2x + 1)— + k'^x^ 

dx^ dx 

+2iskx^ - {I - s){l + s + l)x{x + 1) - isk{2x + 1) + k'^] sRiu> = 0. (El) 

In the region x <^ {I + l)/k (which includes the horizon) the third and fourth terms can 
be neglected, so that the equation has three regular singular points and its general solution 
is expressible in terms of hypergeometric functions: 

sRiu.= Ci 2^1 (-/ - s, / - s + 1; 1 - s - 2tk; -x) x'^-'^x + ly^"" 

+C2 2F1 (-/ + 2ik, 1 + 1 + 2ik- 1 + S + 2ik; -x) {-iyx'''{x + 1)"^+^*= (E2) 

where Ci and C2 are constants. In the region x ^ k + 1 (which stretches to infinity) the 
last two terms of ([ET|) are ignorable and the following general solution of the approximate 
equation can be deduced: 

^Ri^= Di iFi (/ + 1 - s, 2/ + 2; 2ikx) e-^'^V"" 

+D2 iFi (-/ - s, -21; 2ikx) e-'^^^a;-'-"-^ . (E3) 

One can now match these general solutions in the region of overlap /c + 1 -C a; <^ {I + l)/k. 
(The hypergeometric functions in ( |E2|) may be approximated when x k + 1 using Eq. 2 
on p. 108 of Ref. [^, and when x ^ [l + l)/k those in (|E3| ) can be simply replaced by 1.) 
The following matching relations between Ci, C2 and Di, D2 obtain: 
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^ r(2/ + l) T{l-s-m) ^ , , r(2/ + l) T{l + s + 2zk) ^ 

TxT^ITT) TxT+F^ ' ^ r(/ + . + i) r(i + / + 2^/t) ^^^^^ 

r(-2/-i)r(i-.-2^fc) r(-2/-i)r(i + . + 2^A:) 

rn-s) rH-2.A;) ^^ + (-^) r(.-/) rH + 2.A;) ' ^^^^^ 

Consider now the particular solution _2-R|Ji(r) of the s = —2 radial equation. First 
comparing the form of this solution as r — » 2M (see (^T7^)) with that of the general solution 
( p^ ) in the same limit, one observes that 

Ci = {2MYBZ and C2 = (E5) 

for this particular choice. Secondly comparing the form of -2-R;Ji('") with that of ( [E3D as 
r — > cxD (see (|3.7a|) and Eq. 2, p. 278 of Ref. pGf for the asymptotic expansion of the confluent 
hypergeometric function iFi) yields 

^(-2*)-"B. + F^(-«*0-B.^ ^ (E6a) 
^|^(2«)-'+'Z), + St3_(2,k)'»o,= (2A0Mi° . (E6b) 

Eliminating Ci,C2,-Di and 1^2 between (|E4| ) (with s = —2), (|E5|) and (|E6| ) results in a 
pair of equations which determine the incoming reflection and transmission amplitudes; the 
leading behavior of i?;™ as 2Muj — > is then found to be 

A similar comparison of +21^1^ with the approximate hypergeometric solutions yields 

AZ ^ 2(2M)^|^(-4^M.;) (E8) 

for small 2Muo . 



32 



REFERENCES 



* Address until October, 1995. 

t Present address: J.H. Schroder Wagg & Co., 120 Cheapside, London EC2V 6DS, UK 
[1] B. Allen, J. G. Mc Laughlin and A. C. Ottewill, Phys. Rev. D45, 4486 (1992). 
[2] T. Regge and J. A. Wheeler, Phys. Rev. 108, 1063 (1957). 
[3] F. J. Zerilli, Phys.Rev D2, 2141 (1970). 
[4] S. A. Teukolsky, Phys. Rev. Lett. 29, 1114 (1972). 
[5] S. A. Teukolsky, Astrophys. J. 185, 635 (1973). 
[6] P. L. Chrzanowski, Phys. Rev. Dll, 2042 (1975). 

[7] J. B. Hartle and S. W. Hawking, Commun. Math. Phys. 27, 283 (1972). 

[8] P. Candelas, P. Chrzanowski and K. W. Howard, Phys. Rev. D24, 297 (1981). Note 

that our conventions set G = 1 while these authors choose 167rG = 1. 
[9] D. N. Page, Phys. Rev. D13, 198 (1976). 
[10] C. W. Misner, K. S. Thorne and J. A. Wheeler, Gravitation (Freeman, San Francisco, 

1973). 

[11] E. Newman and R. P. Penrose, J. Math. Phys. 3, 566 (1962); E. Newman and R. P. 

Penrose, J. Math. Phys. 7, 863 (1966). Note that our sign conventions are not the same 

as those for these papers — c.f. Appendix A. 
[12] J. Hartle and S. W. Hawking, Commun. Math. Phys. 27, 283 (1973). 
[13] N. H. Barth and S. M. Christensen, Phys. Rev. D28, 1876 (1983). 
[14] R. M. Wald, Phys. Rev. Lett. 41, 203 (1978). 
[15] W. H. Press and S. A. Teukolsky, Astrophys. J., 193, 443 (1974). 
[16] S. M. Christensen and S.A. Fulling, Phys. Rev. D15, 2088 (1977). 
[17] P. Candelas, Phys. Rev. D21, 2185 (1980). 

[18] B. P. Jensen, J. G. Mc Laughlin and A. C. Ottewill, Phys. Rev. D45, 3002 (1992). 
[19] E. Leaver, J. Math. Phys. 27, 1238 (1986). 

[20] S. M. Churilov and A. A. Starobinskii , Zh. Eksp. Teor. Fiz. 65, 3 (1973) [Sov. Phys. 
JETP 38, 1 (1974)]. 

[21] The fact that some of the functions involved in these expansions have a branch point 
at = throws doubt on the validity of this naive analysis. However, a more careful 
analysis based on Appendix A of shows that the result is justified. 

[22] B. Allen, A. Folacci and A. C. Ottewill, Phys. Rev. D38, 1069 (1988). 

[23] S. Chandrasekhar, The Mathematical Theory of Black Holes (Clarendon Press, Oxford, 
1983). 

[24] J. N. Goldberg, A. J. MacFarlane, E. T. Newman, F. Rohrlich and E. C. G. Sudarshan, 
J. Math. Phys. 8, 2155 (1967). 

[25] This formula is essentially a consequence of the group multiplication law. See for exam- 
ple, N.Ja. Vilenkin, Special Functions and the Theory of Group Representation, (Amer- 
ican Mathematical Society, Providence, 1968), p. 131, Eq.(8), with 6*2 = —Oi. We thank 
E. Kalnins for directing us to this reference. 

[26] A. Erdeyli, W. Magnus, F. Oberhettinger, and F. G. Tricomi, Higher Transcendental 
Functions (McGraw-Hill, New York, 1953), Vol. 1. 

[27] I.S. Gradsteyn and LM. Ryzhik, Tables of Integrals, Series and Products, (Academic 
Press, New York, 1980) 



33 



FIGURES 





FIG. 1. Representation of the 'in' and 'up' modes on a Penrose diagram of the exterior 
Schwarzschild space-time. The 'in' modes emanate from T~ and are partial transmitted across 
and partly reflected out to X+. The 'out' modes emanate from H" and are partial transmitted 
out to X"*" and partly reflected across 7^+. 
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FIG. 2. The space-time region V of the exterior Schwarzschild space-time. 






2M 3M 4M 5M^ 

FIG. 7. (2M)6(1 - 2M/r)4(|^'o|^)^~'^ x 10^ The value of this combination at r = 2M is 19-6 
and its asymptotic value at infinity is 15-8. 
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2M 3M 4M 5M^ 

FIG. 8. {2Mf{l-2M/r)-'^{\^i\^)^-^ x 10^ The value of this combination at r = 2M is 9-14 
and its asymptotic value at infinity is 0-988. 



75 



14 



13 



2M 3M 4M 5M^ 

FIG. 9. {2Mf{r/2My^{l-2M/rf{\-^of)^-^ x 10^ The value of this combination at r = 2M 
is 12-1 and its asymptotic value at infinity is 29.1. 
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2M 3M 4M 5M 

FIG. 10. (2M)6(r/2M)2(l-2M/r)(|^4|^)'^--^ X 10"^. The value of this combination at r = 2M 
is 7-54 and its asymptotic value at infinity is 0-0842. 




2M 3M 4M 5M 

FIG. 11. (2M)6(1 - 2M/rf{\i!o\^)^-^ x 10^ = (2M)6l6(l - 2M/r)(|*4|^)^~^ x 10^. The 
value of this combination at r = 2M is 12-1 and its asymptotic value at infinity is 0-158. 
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